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Abstract 


This  report  presents  a  recent  project  for  further  verification  and  enhancement  of  a  recently  developed 
NURBS-based  computational  method  in  VAST.  This  method  has  the  potential  for  close  integration  of 
engineering  analysis  capability  with  Computer  Aided  Design  (CAD)  systems.  It  also  has  the 
advantages  of  preserving  exact  geometry  of  structures  to  be  analyzed  and  providing  a  convenient 
means  for  mesh  refinement.  In  a  previous  project,  the  NURBS-based  formulation  was  implemented 
into  the  VAST  program  and  verified  for  a  wide  variety  of  engineering  analyses,  including  linear  static, 
eigenvalue  and  nonlinear  analyses  involving  plastic  post-collapse  responses  of  thin  shell  structures.  In 
this  work,  additional  verification  of  the  NURBS-based  element  is  performed  and  the  numerical 
problems  observed  in  the  previous  study  were  also  investigated  and  resolved.  A  post-processing 
capability  for  the  NURBS-based  element  is  developed  and  the  possibility  for  applying  this  element  for 
solving  practical  engineering  problems  is  also  explored.  This  study  confirmed  advantages  of  NURBS- 
based  elements  from  the  computational  point  of  view.  An  automatic  model  generator  for  NURBS- 
based  elements  still  needs  to  be  developed  in  the  future. 


Resume 


Ce  rapport  presente  un  projet  recent  de  verification  additionnelle  et  d’ amelioration  d’une  methode  de 
calcul  basee  sur  NURBS  qui  est  integree  a  VAST  et  qui  pourrait  potentiellement  permettre  une 
integration  plus  etroite  des  capacites  d’analyse  d’ingenierie  a  des  systemes  de  conception  assistee  par 
ordinateur  (CAO).  Cette  methode  presente  aussi  l’avantage  de  preserver  la  geometrie  exacte  des 
structures  analysees  et  elle  constitue  un  moyen  pratique  de  raffinement  des  maillages.  Dans  le  cadre 
d’un  projet  precedent,  la  formulation  fondee  sur  NURBS  avait  ete  mise  en  oeuvre  dans  le  programme 
VAST  et  verifiee  avec  une  vaste  ganime  d’analyses  en  ingenierie,  notamment  des  analyses  statiques 
lineaires,  de  valeurs  propres  et  non  lineaires  portant  sur  les  reactions  plastiques  apres  effondrement  de 
structures  a  enveloppe  mince.  Au  cours  de  ces  travaux,  une  verification  plus  poussee  de  Telement  base 
sur  NURBS  a  ete  effectuee  et  les  problemes  numeriques  observes  au  cours  de  T  etude  precedente  ont 
aussi  ete  etudies  et  resolus.  Une  capacite  de  post-traitement  est  en  cours  de  developpement  et  on 
explore  la  possibilite  d’appliquer  cet  element  a  la  resolution  de  problemes  d’ingenierie  concrets.  Cette 
etude  a  confirme  les  avantages  de  Telement  base  sur  NURBS  du  point  de  vue  computationnel.  Par 
ailleurs,  il  est  necessaire  de  developper  un  generateur  de  modeles  automatise  pour  des  elements  bases 
sur  NURBS. 
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Executive  summary 


Development  of  a  NURBS-based  Computational  Method  for  the 
VAST  Finite  Element  Program:  Phase  2  -  Final  Report 

Lei  Jiang,  Eric  Teng;  DRDC  Atlantic  CR  2008-118;  Defence  R&D  Canada  -  Atlantic; 
October  2008 

Introduction:  A  new  numerical  method  for  solving  problems  governed  by  partial  differential 
equations,  such  as  solids  and  structures,  has  been  developed  recently,  taking  inspiration  from  Computer 
Aided  Design  (CAD).  Although  this  method  has  many  features  in  common  with  the  finite  element 
method  and  some  features  in  common  with  the  meshless  methods,  it  is  more  geometrically  based  and 
provides  a  potential  for  a  close  integration  with  CAD  systems.  This  numerical  formulation  is  based  on 
NURBS  (Non-Uniform  Rational  B-Splines),  a  standard  technology  employed  in  CAD  systems.  Using 
this  method,  exact  geometric  representations  of  the  solids  and  structures  are  maintained  for  any  level  of 
mesh  refinement. 

Results:  The  NURBS-based  computational  method  was  implemented  into  the  VAST  finite  element 
program  in  a  previous  contract  as  a  special  element  type,  similar  to  the  conventional  solid  element. 
Under  the  present  contract,  further  verification  of  the  NURBS-based  element  was  performed  and  the 
performance  of  the  NURBS-based  and  traditional  20-noded  solid  elements  was  extensively  compared. 
The  previously  observed  convergence  difficulty  of  the  NURBS-based  element  was  investigated  and 
resolved.  To  improve  user  friendliness  of  the  NURBS-based  element,  a  post-processing  capability  was 
developed.  Finally,  a  NURBS-based  finite  element  model  of  a  simplified  ship  hull  structural 
component  consisting  of  a  hull  plate  and  intersecting  truncated  transverse  bulkhead  and  two 
longitudinal  stiffeners  was  created  and  analyzed.  The  results  indicated  that  to  achieve  the  same  level  of 
accuracy,  the  size  of  the  NURBS-based  finite  element  model  was  considerably  smaller  than  that  of  the 
traditional  20-noded  solid  element  model,  resulting  in  a  much  shorter  computation  time. 

Significance:  The  present  study  confirmed  the  advantages  of  the  NURBS-based  element  formulation 
from  the  computational  point  of  view.  This  formulation  provides  an  exact  representation  of  the  problem 
geometry  and  a  convenient  approach  for  further  mesh  refinement  through  insertion  of  knots  and  control 
points.  The  NURBS-based  element  is  now  working  reliably  for  all  types  of  analyses  and  its  potential 
for  analyses  of  ship  structural  components  has  been  demonstrated. 

Future  plans:  Significant  developments  are  still  required  before  it  can  be  utilized  in  practical 
engineering  analyses  of  complete  ship  structures  and  structural  connection  details.  A  pre-processing 
capability  is  needed  to  take  the  CAD  model  and  transform  it  into  an  initial  mesh  or  control  net  of 
NURBS-based  solid  finite  elements.  In  addition,  the  possibility  for  developing  NURBS-based  shell 
elements  should  be  explored. 
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Introduction  :  Une  nouvelle  methode  numerique  de  resolution  de  problemes  fondee  sur  des  equations 
differentielles  partielles,  comme  celles  appliquees  aux  solides  et  aux  structures,  a  ete  developpee 
recemment  en  s’ inspirant  de  la  conception  assistee  par  ordinateur  (CAO).  Bien  que  cette  methode 
partage  plusieurs  caracteristiques  avec  la  methode  utilisant  des  elements  finis  et  quelques 
caracteristiques  avec  les  methodes  sans  maillage,  elle  est  plutot  basee  sur  la  geometrie  et  elle  presente 
le  potentiel  d’une  integration  etroite  avec  les  systemes  de  CAO.  Cette  formulation  mathematique  est 
basee  sur  NURBS  (courbes  B-splines  rationnelles  non  uniformes),  une  technologie  standard  utilisee 
dans  les  systemes  de  CAO.  Au  moyen  de  ce  systeme,  il  est  possible  de  conserver  des  representations 
geometriques  exactes  de  solides  et  de  structures  pour  tout  degre  de  raffmement  des  maillages  et  aucune 
communication  ulterieure  avec  la  description  CAO  n’est  requise  au  cours  du  processus  de  raffmement 
du  maillage  au  moyen  d’elements  finis. 

Resultats  :  La  methode  de  calcul  basee  sur  NURBS  a  ete  mise  en  oeuvre  dans  le  programme  de  calcul 
par  elements  finis  VAST  dans  le  cadre  d’un  contrat  precedent  sous  forme  d’un  type  d’element  special, 
semblable  a  un  element  solide  classique.  Dans  le  cadre  du  present  contrat,  on  a  effectue  des 
verifications  plus  poussees  de  Telement  base  sur  NURBS  et  on  a  effectue  une  comparaison  poussee  des 
elements  solides  bases  sur  NURBS  et  des  solides  classiques  a  20  noeuds.  Le  probleme  de  convergence 
touchant  Telement  base  sur  NURBS,  releve  precedemment,  a  ete  etudie  et  resolu.  Afm  d’ameliorer  la 
convivialite  de  Telement  base  sur  NURBS,  une  capacite  de  post-traitement  a  ete  developpee. 
Finalement,  on  a  cree  et  analyse  un  modele  d’element  fmi  base  sur  NURBS  representant  une 
composante  structural  de  coque  de  navire  simplifiee  composee  d’un  borde  de  carene  (tole  de  coque), 
d’une  cloison  transversale  tronquee  en  intersection  et  de  deux  lisses  longitudinales,  a  ete  creee  et 
analysee.  Les  resultats  obtenus  indiquent  que,  pour  obtenir  le  meme  degre  de  precision,  la  taille  du 
modele  a  elements  finis  base  sur  NURBS  etait  considerablement  plus  petite  que  celle  du  modele  a 
element  solide  a  20  noeuds  classique,  ce  qui  reduit  considerablement  le  temps  de  calcul. 

Importance  :  L’etude  actuelle  a  confirme  les  avantages  de  la  formulation  d'elements  basee  sur  NURBS 
du  point  de  vue  de  la  quantite  de  calcul  requise.  Cette  formulation  foumit  une  representation  exacte  de 
la  geometrie  du  probleme  et  elle  constitue  une  approche  commode  pour  raffiner  le  maillage  en  inserant 
des  noeuds  et  des  points  de  controle.  L’element  base  sur  NURBS  fonctionne  de  maniere  Liable  pour 
tous  les  types  d’analyse  et  nous  avons  demontre  son  potentiel  dans  le  cas  de  Tanalyse  des  composantes 
structurelles  des  navires. 

Plans  futurs  :  Des  travaux  de  developpement  additionnels  sont  encore  requis  afin  que  cette  methode 
puisse  etre  appliquee  a  des  analyses  d’ingenierie  pratiques  de  structures  completes  de  navires  et 
d’ analyses  des  details  des  assemblages  structuraux.  Une  fonction  de  pretraitement  est  requise  afm  de 
pouvoir  transformer  les  modeles  de  CAO  en  des  maillages  initiaux  ou  des  reseaux  de  controle 
d’elements  finis  solides  bases  sur  NURBS.  En  outre,  il  serait  necessaire  d’explorer  la  possibilite  de 
developper  des  elements  d’enveloppe  bases  sur  NURBS. 


tv 
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INTRODUCTION 


One  of  the  goals  of  the  research  project  entitled  “Investigation  of  LCM  through  Single  Product 
Models  (SPM)”  is  to  automate  the  process  of  generation  of  structural  finite  element  models 
from  geometric  information  in  SPM  databases.  NURBS  (Non-Uniform  Rational  B-Spline)  are 
a  standard  technology  employed  in  SPM/CAD  systems  to  represent  geometry  of  surfaces  such 
as  ship  hull  forms.  A  numerical  method  for  solving  problems  governed  by  partial  differential 
equations,  such  as  stress  analyses  of  solids  and  structures,  was  developed  recently  based  on 
NURBS  representation  of  structural  components.  Although  this  method  has  many  features  in 
common  with  the  traditional  finite  element  method,  and  some  features  in  common  with  the 
meshless  methods,  it  is  more  geometrically  based  and  has  potential  to  directly  utilize  NURBS- 
based  geometry  definition  in  SPM/CAD  systems  to  perfonn  engineering  analyses.  It  thus 
offers  a  number  of  advantages  over  the  traditional  finite  element  methods  for  integration  with 
SPM  databases. 

The  initial  phase  of  this  work  involved  implementation  and  verification  of  this  NURBS-based 
computational  method  in  the  VAST  finite  element  program  [1].  The  implementation  was 
achieved  as  a  new  element  type,  IEC82,  where  the  input  data  included  the  polynomial  orders, 
the  knots  vectors,  the  control  points  and  the  weights.  The  verification  problems  included  linear 
static,  natural  frequency,  linear  buckling  and  nonlinear  analyses  involving  both  geometric  and 
material  nonlinearities.  The  numerical  results  demonstrated  a  number  of  advantages  of  this 
NURBS-based  formulation  from  the  computational  point  of  view.  However,  the  test  problems 
considered  previously  in  [1]  only  involved  relatively  simple  geometries,  such  as  cylindrical 
and  spherical  shells,  and  no  stiffeners  were  considered  [1].  In  addition,  the  convergence  study 
performed  in  the  previous  phase  was  incomplete  because  the  reference  solutions  employed  in 
the  study  were  obtained  using  a  relatively  coarse  mesh  of  20-noded  solid  elements.  For  linear 
buckling  analyses,  significant  discrepancies  were  noticed  between  buckling  loads  and  mode 
shapes  predicted  by  the  NURBS-based  and  conventional  solid  elements.  Finally,  the  efficiency 
of  the  NURBS-based  method  needs  to  be  tested  using  a  more  complex  ‘ship-like’  structural 
component  and  a  post-processing  capability  is  still  to  be  developed  for  viewing  model  meshes 
and  results. 

The  present  phase  of  the  project  is  to  address  these  remaining  issues  from  the  previous  phase 
of  work,  as  outlined  above.  The  following  work  items  are  to  be  perfonned: 

1.  Conduct  a  more  complete  convergence  study  of  the  NURBS-based  finite  element  for 
linear  static  analysis  and  compare  the  convergence  properties  of  the  NURBS-based 
solid  element  with  that  of  the  traditional  20-noded  solid  element. 

2.  Investigate  the  accuracy  issue  of  the  NURBS-based  solid  element  for  linearized 
buckling  analysis,  through  a  detailed  review  of  the  mathematical  formulation  and 
computer  code  and  additional  numerical  verifications. 

3.  Compare  the  computational  efficiency  of  the  traditional  20-noded  solid  and  NURBS- 
based  solid  elements  for  practical  engineering  problems.  This  task  includes  generation 
of  20-noded  and  NURBS-based  solid  element  models  of  simplified  ship  structural 
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components  and  detailed  comparison  of  efficiency  and  accuracy  of  these  models  for 
linear  static  analysis. 

4.  Develop  a  post-processing  capability  to  automatically  generate  dummy  20-noded  solid 
element  models  in  knots  space  and  to  map  and  display  displacement  and  stress  results 
on  to  the  dummy  solid  element  models,  so  these  20-noded  solid  element  models,  along 
with  the  displacements  and  stress  results,  can  be  readily  displayed  with  the  existing 
graphics  capability  provided  in  Trident  FEA. 

All  these  tasks  have  been  accomplished  successfully  and  the  details  will  be  presented  in  the 
following  chapters  in  this  report.  The  conclusions  from  this  study  are  given  in  Chapter  6. 
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2.0  DEVELOPMENT  OF  A  POST-PROCESSING  CAPABILITY  FOR  THE 
NURBS-BASED  SOLID  ELEMENT 

The  development  and  verification  of  the  post-processing  capability  is  presented  in  this  chapter. 
Although  this  is  the  last  work  item  in  the  Statement  of  Work,  we  decided  to  discuss  it  first  in 
the  report  because  this  capability  is  extensively  utilized  in  the  study  of  other  work  items  in  this 
contract. 

As  presented  in  the  previous  contract  report  [1],  the  “element”  in  the  NURBS-based  solid 
element  formulation  is  defined  by  knot  intervals  in  the  knots  space.  However,  differing  from 
the  traditional  finite  elements,  the  approximation  space  for  the  NURBS-based  element  is  non¬ 
local,  which  means  that  the  displacement  field  in  an  element  depends  not  only  on  degrees-of- 
freedom  of  nodes  inside  the  domain  of  the  element,  but  also  on  nodes  outside  the  element. 
This  feature  is  in  common  with  the  meshless  or  mesh-free  method. 

In  the  post-processing  capability,  each  “element”  in  a  NURBS-based  solid  element  model  is 
subdivided  into  a  uniform  mesh  of  20-noded  solid  elements.  This  subdivision  is  first 
performed  in  the  knot  space  and  then  the  actual  coordinates  of  the  nodes  of  the  20-noded  solid 
elements  are  calculated  from  the  coordinates  of  the  control  nodes  in  the  three-dimensional 
physical  space  using  NURBS  functions.  The  numbers  of  elements  in  all  three  directions,  used 
for  this  subdivision,  are  user  controlled  parameters,  so  a  finer  subdivision  can  always  be 
conveniently  achieved  when  higher  resolutions  are  desired  for  result  presentation. 

When  all  “elements”  in  the  NURBS-based  model  are  subdivided,  node  equivalencing  is 
required  to  eliminate  the  redundant  nodes  at  the  element  interfaces.  This  operation  is  time- 
consuming,  but  is  required  to  ensure  stress  contour  plots  can  be  displayed  properly.  At  the 
present  time,  the  post-processing  capability  is  able  to  process  results  from  both  linear  and 
nonlinear  analyses,  including  displacements,  stresses  and  plastic  strains.  Among  them,  the 
displacements  are  calculated  using  the  displacements  at  the  control  points  through  the  NURBS 
functions,  just  like  the  nodal  coordinates. 

The  calculation  of  stresses  and  strains  in  the  dummy  20-noded  solid  element  model  is  a  little 
different.  In  VAST,  element  stresses  in  the  20-noded  solid  element  are  stored  at  the  8  corner 
nodes  and  this  fonnat  of  stress  storage  is  also  adopted  in  the  NURBS-based  solid  element. 
During  post-processing,  the  stress  variation  in  a  NURBS-based  element  is  approximated  by  a 
tri-linear  function  from  which  the  stresses  at  comers  of  each  element  in  the  dummy  20-noded 
element  model  can  be  computed.  In  nonlinear  analyses  involving  plastic  deformations,  stresses 
and  plastic  strains  are  calculated  at  the  integration  points.  In  this  situation,  the  eight  outmost 
integration  points  are  used  for  the  tri-linear  interpolation.  Alternative  interpolation  procedures, 
including  higher-order  interpolation  functions  and  piecewise  linear  interpolations,  have  been 
explored.  However,  the  results  are  not  satisfactory. 

The  post-processing  program  is  developed  in  FORTRAN  in  a  style  consistent  with  the  VAST 
code.  As  a  result,  it  can  be  conveniently  incorporated  into  VAST  when  it  is  required  in  the 
future. 
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In  order  to  demonstrate  and  verify  this  post-processing  capability,  two  example  problems 
considered  in  the  previous  phase  [1]  were  utilized.  These  test  problems  involved  geometrically 
nonlinear  analysis  of  a  hinged  cylindrical  shell  and  combined  geometrically  and  materially 
nonlinear  analysis  of  a  hinged  spherical  shell,  both  subjected  to  concentrated  loads. 

The  control  net  of  the  NURBS-based  solid  element  model,  the  dummy  20-noded  solid  element 
models  created  by  the  post-processor  and  the  deformation  contour  plot  generated  by  Trident 
using  the  dummy  20-noded  solid  element  model  are  shown  in  Figures  2. 1-2.3  and  2.4-2. 6, 
respectively  for  these  two  test  cases.  These  figures  indicate  that  the  nodal  displacements  of  the 
dummy  20-noded  solid  element  model  are  computed  correctly.  It  is  also  interesting  to  note  that 
although  the  outer  edge  of  the  control  net  for  the  second  test  problem  is  non-smooth,  the  20- 
noded  solid  element  model  always  represent  the  exact  geometry  of  the  structure.  Actually,  this 
is  one  of  the  advantages  for  using  NURBS  representation  to  describe  a  problem  geometry. 

For  verification  of  the  post-processing  capability  for  stresses  and  strains,  we  utilized  the 
second  test  case  mentioned  above,  which  involved  a  combined  geometrically  and  materially 
nonlinear  analysis  of  a  hinged  spherical  shell  subjected  to  a  centre  point  load.  As  shown  in 
Figure  2.4,  in  the  NURBS-based  solid  element  model,  one  quarter  of  the  spherical  shell  was 
represented  by  three  NURBS-patches.  Each  of  these  patches  contains  a  4x4x1  mesh  of 
elements  defined  on  knots  intervals  and  during  the  post-processing,  each  element  is  subdivided 
into  a  4x4x10  mesh  of  20-noded  solid.  The  reason  for  using  ten  elements  in  the  thickness 
direction  is  to  adequately  capture  through-the-thickness  variations  of  stresses  and  plastic 
strains  using  Trident  FEA,  where  the  stress  and  strain  are  displayed  based  on  values  at  element 
centroid.  The  generated  stress  and  plastic  strain  contour  plots  for  a  typical  load  step  are  given 
in  Figures  2.7  and  2.8,  in  which  through-the-thickness  stress  and  strain  gradients  are  observed. 
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Figure  2.1:  Control  Net  of  NURBS-based  Solid  Element  Model  of  a  Cylindrical  Shell 


Figure  2.2:  A  20-Noded  Solid  Element  Model  of  the  Cylindrical  Shell  Generated  by  the 

Post-Processor 
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Figure  2.3:  Displacement  Contour  Generated  Using  the  Post-Processing  Capability  for 
Geometrically  Nonlinear  Analysis  of  a  Hinged  Cylindrical  Shell 
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Figure  2.4:  Control  Net  of  NURBS-based  Solid  Element  Model  of  a  Spherical  Shell 


Figure  2.5:  A  20-Noded  Solid  Element  Model  of  the  Spherical  Shell  Generated  by  the 

Post-Processor 


DRDC  Atlantic  CR  2008-118 


7 


Displacement  Contours 


Displacements;  Load  Step  50  Load  Parameter  0.491  E+01 


X  Translation 

3.83E-04 
-6.59E-03 
-1.36E-02 
-2.05E-02 
-2.75E-02 
-3.45E-02 
-4.14E-02 
-4.84E-02 
-5.54E-02 
-6.23E-02 
-6.93E-02 
-7.63E-02 
-8.32E-02 
-9.02E-02 
-9.72E-02 
-1.04E-01 
-1.1  IE-01 
-1.18E-01 
-1.25E-01 


Figure  2.6:  Displacement  Contour  Generated  Using  the  Post-Processing  Capability  for 
Geometrically  Nonlinear  Analysis  of  a  Hinged  Spherical  Shell 


8 


DRDC  Atlantic  CR  2008-118 


Contour  Plot 
Stress  Components 
Load  Step  8 
20N/Solid  SYY 


2.01  E+04 
1.77E+04 
1.53E+04 
1.29E+04 
1.05E+04 
8.13E+03 
5.73E+03 
3.33E+03 
9.35E+02 
-1.46E+03 
-3.86E+03 
-6.26E+03 
-8.66E+03 
-1.11  E+04 
-1.35E+04 
-1.59E+04 
-1.82E+04 
-2.06E+04 
-2.30E+04 


Contour  Plot 
Stress  Components 
Load  Step  8 
20N/Solid  SYY 


2.01  E+04 
1.77E+04 
1.53E+04 
1.29E+04 
1.05E+04 
8.13E+03 
5.73E+03 
3.33E+03 
9.35E+02 
-1.46E+03 
-3.86E+03 
-6.26E+03 
-8.66E+03 
-1 .1 1  E+04 
-1.35E+04 
-1.59E+04 
-1.82E+04 
-2.06E+04 
-2.30E+04 


Figure  2.7:  Stress  Contour  Plots  Generated  Using  the  Post-Processor  at  a  Typical  Load 
Step  in  Nonlinear  Analysis  of  a  Hinged  Spherical  Shell 
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Figure  2.8:  Plastic  Strain  Contour  Plots  Generated  Using  the  Post-Processor  at  a  Typical 
Load  Step  in  Nonlinear  Analysis  of  a  Hinged  Spherical  Shell 
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3.0  CONVERGENCE  STUDY  OF  THE  NURBS-BASED  SOLID  ELEMENT  FOR 
LINEAR  STATIC  ANALYSIS 

In  the  previous  phase  of  development  of  the  NURBS-based  solid  element  for  VAST  [1],  the 
convergence  characteristics  of  this  element  formulation  were  studied  by  considering  linear 
elastic  response  of  a  hinged  cylindrical  shell  subjected  to  a  centre  point  load.  Due  to  symmetry 
in  the  problem,  only  one  quarter  of  the  shell  needed  to  be  considered  in  the  numerical  solution. 
In  the  previous  work  [1],  a  variety  of  NURBS-based  solid  element  models  were  constructed 
using  different  polynomial  orders  and  different  numbers  of  knot  intervals.  The  displacement 
and  stress  results  obtained  by  the  NURBS-based  solid  element  were  compared  against 
solutions  of  20-noded  solid  element  based  on  the  conventional  isoparametric  fonnulation. 
However,  at  that  time,  only  one  relatively  coarse  mesh  of  the  20-noded  solid  element  was 
considered,  so  the  convergence  property  of  the  NURBS-based  element  was  not  fully  verified. 

In  the  present  work,  the  previous  convergence  study  was  extended  by  refining  the  reference 
solutions  using  additional  meshes  of  the  conventional  20-noded  solid  element.  These  20-noded 
solid  element  meshes  could  be  divided  into  three  groups  based  on  their  general  layout.  The 
meshes  in  the  first  group  could  be  represented  as  NxNx2  and  were  used  to  evaluate  the  effect 
of  mesh  refinement  in  the  in-plane  directions.  The  second  group  contained  meshes  defined  as 
4x4xN  and  were  used  to  reveal  the  effect  of  mesh  refinement  through-the-thickness  direction. 
The  final  group  contained  a  number  of  highly  refined  meshes,  which  were  utilized  to  provide 
references  for  comparison.  Four  typical  meshes  of  the  20-noded  solid  element  used  in  the 
convergence  study  are  displayed  in  Figure  3.1  and  the  results  are  summarized  in  Table  3.1  for 
both  full  and  reduced  integrations. 

The  newly  obtained  results  for  the  20-noded  solid  element  were  presented  in  Table  3.1  and 
also  compared  with  the  solutions  of  the  NURBS-based  solid  element  in  Figures  3.2  and  3.3, 
where  the  results  of  the  NURBS-based  solid  elements  are  taken  from  the  previous  report  [1]. 
These  comparisons  indicate  that  for  models  of  the  same  number  of  nodes,  the  NURBS-based 
solid  element  is  far  more  accurate  than  the  conventional  20-noded  solid  element.  For  example, 
the  finest  NURBS-based  element  model  considered  in  the  previous  study  contained 
13x13x2=338  control  points  and  predicted  a  centre  deflection  of  2.785.  On  the  other  hand,  the 
centre  deflection  predicted  by  a  5x5x2  mesh  (360  nodes)  of  the  20-noded  solid  element  is  only 
2.487  for  full  integration  and  2.597  when  reduced  integration  was  utilized.  The  prediction  of 
the  NURBS-based  solid  element  was  in  a  far  better  agreement  with  the  reference  solution  of 
2.848  obtained  using  a  30x30x10  mesh  (40641  nodes)  of  the  20-noded  solid  element. 

In  addition,  these  results  also  indicated  that  for  this  particular  test  case,  the  accuracy  of  the  20- 
noded  solid  element  is  very  sensitive  to  the  order  of  numerical  integration.  Use  of  full 
integration  results  in  “locking”  of  the  element,  so  the  convergence  rate  becomes  extremely 
slow.  The  results  of  the  NURBS-based  element  seem  to  be  less  sensitive  to  the  numerical 
integration  orders.  All  results  for  the  NURBS-based  element  presented  in  Figures  3.2  and  3.3 
were  obtained  using  full  integration. 
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Table  3.1:  Results  of  Convergence  Study  Using  20-Noded  Solid  Element  -  Centre 
Deflection  of  a  Hinged  Cylindrical  Shell  Subjected  to  a  Point  Load 


Mesh 

No.  of  Nodes 

Full  Integration 

Reduced  Integration 

Solid  Element  Mesh  NxNx2 

3x3x2 

152 

2.314 

2.606 

4x4x2 

245 

2.430 

2.596 

5x5x2 

360 

2.487 

2.597 

6x6x2 

497 

2.520 

2.602 

7x7x2 

656 

2.541 

2.608 

8x8x2 

837 

2.556 

2.614 

9x9x2 

1040 

2.568 

2.620 

Solid  Element  Mesh  4x4xN 

4x4x2 

245 

2.430 

2.565 

4x4x4 

425 

2.434 

2.639 

4x4x6 

605 

2.435 

2.699 

4x4x8 

785 

2.435 

2.765 

4x4x10 

96 

2.435 

2.836 

Refined  Solid  Element  Meshes 

9x9x8 

3320 

2.577 

2.697 

15x15x10 

10656 

2.628 

2.752 

30x30x10 

40641 

2.717 

2.848 

12 
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Figure  3.1:  Typical  20-Noded  Solid  Element  Meshes  Utilized  in  Convergence  Study 
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Center  Deflection  17  Center  Deflection 


Comparison  of  Centre  Deflections  Predicted  by  NURBS-based  Solid  Element 
and  20-Noded  Solid  Element  with  Full  Numerical  Integration 


Figure  3.3:  Comparison  of  Centre  Deflections  Predicted  by  NURBS-based  Solid  Element 
and  20-Noded  Solid  Element  with  Reduced  Numerical  Integration 
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AN  INVESTIGATION  OF  ACCURACY  OF  THE  NURBS-BASED  SOLID 
ELEMENT  FOR  BUCKLING  ANALYSIS 


In  the  previous  phase  of  the  development  of  NURBS-based  computational  method  for  VAST 
[1],  a  test  case  involving  linear  buckling  analyses  of  a  hinged  cylindrical  shell  subjected  to  a 
centre  point  load  was  considered.  The  problem  geometry  and  loading  condition  are  shown  in 
Figure  4.1,  where  hinged  boundary  conditions  are  applied  along  the  straight  edges  and  the 
curved  edges  are  stress  free.  Three  NURBS-based  models,  which  all  contain  the  same  number 
of  elements  but  different  orders  of  polynomial  functions,  were  utilized  in  the  verification.  The 
control  net  for  one  of  the  NURBS-based  solid  element  models  is  depicted  in  Figure  4.2,  on 
which  the  boundary  conditions  and  external  load  are  also  shown.  During  the  earlier  analyses, 
the  minimum  orders  of  numerical  integration,  which  produced  non-singular  stiffness  matrices, 
were  employed  and  the  predicted  buckling  load  factors  for  the  first  5  modes  were  compared  to 
the  results  of  a  20-noded  solid  element  model.  The  figure  showing  this  comparison  is  copied 
from  the  previous  report  [1]  and  presented  in  Figure  4.3.  These  results  indicated  that  the 
second  and  fifth  modes  obtained  by  the  NURBS-based  solid  elements  are  not  converging  with 
the  increase  of  the  polynomial  order. 

In  order  to  better  understand  the  problem,  we  first  compared  the  buckling  mode  shapes  for 
Modes  2  and  5  predicted  by  the  various  models.  The  buckling  mode  shapes  obtained  using  the 
20-noded  solid  element  and  the  NURBS-based  solid  element  models  with  the  third-  and  fifth- 
order  polynomials  are  given  in  Figures  4.4,  4.5  and  4.6,  respectively.  These  results  indicated 
that  the  mode  shapes  predicted  by  the  20-noded  solid  and  the  NURBS-based  element  with  3ld- 
order  basis  functions  are  consistent,  but  the  buckling  modes  predicted  by  the  NURBS-based 
model  with  higher  order  polynomials  are  very  different. 

To  identify  potential  errors  introduced  during  the  implementation  of  the  NURBS-based  solid 
element,  the  source  code  of  VAST  was  then  checked  carefully,  with  special  attention  paid  to 
the  sections  for  generating  geometric  stiffness  matrix  and  calculating  stresses  at  the  numerical 
integration  points  from  the  comer  values.  However,  no  bugs  were  discovered. 

In  order  to  gain  more  insight  of  this  problem,  it  is  desirable  to  consider  a  simpler  test  example 
having  known  analytical  solutions.  To  this  end,  a  simply-supported  flat  square  plate  subjected 
to  uniform  bi-axial  compression,  as  shown  in  Figure  4.7,  was  utilized.  For  this  problem,  the 
critical  compressive  stresses  at  the  edges  are  expressed  as  [2] 


n  D  [  2  .  2  \ 
cr„  = — —I  m  +  n  I 


a  h 


where  a  and  h  indicate  the  edge  length  and  thickness  of  the  square  plate,  and  D  is  the  flexural 
rigidity  of  plate  defined  as  D=Eh3/\2(\- v^).  Parameters  in  and  n  detennine  the  number  of  half¬ 
waves  in  the  X  and  Y  directions. 


Four  different  models  of  the  NURBS-based  solid  element  were  created  to  represent  one 
quarter  of  the  square  plate.  These  models  contained  the  same  numbers  of  control  points,  but 
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different  polynomial  orders.  As  a  result,  the  number  of  “elements”,  which  were  defined  by 
knot  intervals,  varied  from  model  to  model  because  the  size  of  knot  vectors  was  different  in 
these  models.  The  control  net  of  a  representative  NURBS-based  model  is  also  given  in  Figure 
4.7. 

Numerical  calculations  were  first  performed  using  the  NURBS-based  model  with  the  highest 
polynomial  order  (5x5x2),  but  different  orders  of  numerical  integration.  The  predicted  critical 
stresses  for  the  first  ten  buckling  modes  by  using  different  integration  orders  are  compared  in 
Table  4.1,  where  the  corresponding  analytical  solutions  are  also  shown.  The  numbers  of  half¬ 
waves,  m  and  n,  were  detennined  by  visually  examining  the  buckling  mode  shapes  using  the 
post-processing  capability  presented  in  the  preceding  chapter  of  this  report.  The  mode  shapes 
for  two  representative  buckling  modes  are  shown  in  Figure  4.8,  which  correspond  to  m=n= 3 
and  m=n= 5,  respectively.  The  results  presented  in  Table  4.1  and  Figure  4.8  suggested  that  to 
achieve  reliable  and  accuracy  buckling  loads,  the  orders  of  numerical  integration  should  be  at 
least  identical  to  the  polynomial  orders  used  in  the  NURBS-based  solid  element  model.  These 
results  also  indicated  that  the  numerically  predicted  critical  stresses  are  consistently  lower  than 
the  analytical  solutions.  This  is  because  the  analytical  solutions  are  derived  using  the  plate 
theory  in  which  through-the-thickness  deformations  are  neglected. 

Based  on  the  results  presented  in  Table  4.1,  the  VAST  solver  was  modified  to  set  the  default 
orders  of  numerical  integration  equal  to  the  polynomial  orders  in  each  of  the  parametric 
directions.  To  confirm  the  robustness  of  this  default  integration  rule,  additional  calculations 
were  carried  out  using  other  NURBS-based  models  using  the  default  integration  orders.  The 
numerical  results  are  summarized  in  Table  4.2,  where  a  stable  convergence  pattern  is  observed. 

At  this  point,  the  buckling  problem  for  the  hinged  cylindrical  shell,  depicted  in  Figure  4.1,  was 
reanalyzed  using  the  default  numerical  integration  orders.  The  new  results  of  buckling  load 
factors  for  the  first  five  modes  are  shown  in  Figure  4.9,  where  convergence  is  observed  for  all 
the  modes.  The  buckling  mode  shapes  obtained  by  the  NURBS-based  solid  element  model 
with  5th-order  polynomials  using  the  default  numerical  integration  orders  are  shown  in  Figure 
4.10.  Comparing  these  buckling  mode  shapes  with  those  shown  in  earlier  figures,  Figures  4.4- 
4.6,  it  is  readily  recognized  that  the  errors  observed  in  the  previous  study  [1]  were  caused  by 
the  inadequacy  of  the  numerical  integrations. 

In  order  to  see  whether  this  inadequacy  can  be  observed  by  examining  the  solutions  of  linear 
static  analyses,  we  compared  displacement  and  stress  contours  produced  by  the  same  NURBS- 
based  model,  but  with  different  numerical  integration  orders  in  Figures  4.11-4.12. 
Surprisingly,  these  results  are  in  very  close  agreement.  These  observations  suggest  that  higher 
numerical  integration  orders  are  required  to  accurately  integrate  the  geometric  stiffness  matrix 
than  required  for  the  material  stiffness  matrix. 
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Table  4.1:  Comparison  of  Critical  Stresses  for  Buckling  of  a  Bi- Axially  Loaded  Simply 
Supported  Square  Obtained  by  a  NURBS-based  Solid  Element  Model  Using  Different 

Orders  of  Numerical  Integration 


Mode  # 

NURBS-based  Fini 
Numerical  Integral 

te  Elements  (MPa): 
ion  Order  nuxnvxns 

Analytical  (MPa) 

3x3x2 

4x4x2 

5x5x2 

8x8x3 

m 

n 

Critical  a 

1 

58.37 

58.87 

58.89 

58.89 

1 

1 

59.31 

2 

293.47 

294.41 

294.45 

294.45 

3 

1 

296.56 

3 

293.49 

294.42 

294.45 

294.46 

1 

3 

296.56 

4 

520.94 

525.70 

525.86 

525.88 

3 

3 

533.81 

5 

759.97 

761.11 

761.18 

761.21 

6 

759.99 

761.13 

761.19 

761.23 

7 

979.93 

987.14 

987.39 

987.45 

5 

3 

1008.31 

8 

980.13 

987.26 

987.50 

987.55 

3 

5 

1008.31 

9 

1425.2 

1438.70 

1439.10 

1439.30 

5 

5 

1482.81 

10 

1454.2 

1467.40 

1468.80 

1469.00 

Table  4.2:  Comparison  of  Critical  Stresses  for  Buckling  of  a  Bi- Axially  Loaded  Simply 
Supported  Square  Obtained  by  Different  NURBS-based  Solid  Element  Models  Using  the 

Default  Numerical  Integration 


Mode  # 

NURBS-based  Fini 
Polynomial  O 

te  Elements  (MPa): 

>rder  puxpvxps 

Analytical  (MPa) 

2x2x2 

3x3x2 

4x4x2 

5x5x2 

m 

n 

Critical  a 

1 

59.00 

58.96 

58.91 

58.89 

1 

1 

59.31 

2 

296.24 

294.77 

294.50 

294.45 

3 

1 

296.56 

3 

296.25 

294.78 

294.51 

294.45 

1 

3 

296.56 

4 

528.72 

526.69 

526.10 

525.86 

3 

3 

533.81 

5 

797.99 

768.73 

762.31 

761.18 

6 

798.00 

768.74 

762.32 

761.19 

7 

1017.80 

993.93 

988.52 

987.39 

5 

3 

1008.31 

8 

1017.90 

994.02 

988.63 

987.50 

3 

5 

1008.31 

9 

1481.50 

1448.20 

1440.80 

1439.10 

5 

5 

1482.81 

10 

1738.00 

1545.90 

1487.00 

1468.80 
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R  =  2540  mm 
1  =  254  mm 

h  =  12.7  mm,  6.35  mm 
0  =  0.1  rad.  0 
E  =  3.10275  kN/mm^ 
v  =  0 . 3 


Figure  4.1:  Problem  Definition  for  Buckling  Analyses  of  a  Hinged  Cylindrical  Shell 

Subjected  to  a  Centre  Point  Load 
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Boundary  Conditions 


Force  Resultant;  LoadCase  1: 

STATIC  LOAD  ANALYSIS:  CASE  NUM  BER  1 


Figure  4.2:  Control  Net  of  a  NURBS-based  Solid  Element  Model  Representing  One 
Quarter  of  the  Cylindrical  Shell  with  Boundary  Conditions  and  External  Load 
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Figure  4.3:  Buckling  Loads  Obtained  by  NURBS-based  Solid  Elements  of  Various 
Orders  and  20-Noded  Solid  Elements  Presented  in  Previous  Contract  Reported  [1] 
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Mode  Shape  Contours 


Buckling  Mode  Shape  2  Critical  Load  Factor  0.1468E+02 


X  Translation 

4.1 5E+00 
3.84E+00 
3.53E+00 
3.22E  +  00 
2.90E  +  00 
2.59E+00 
2.28E+00 
1.97E  +  00 
1.66E+00 
1.34E  +  00 
1.03E  +  00 
7.21  E-01 
4.09E-01 
9.77E-02 
-2.14E-01 
-5.26E-01 
-8.38E-01 
-1.15E  +  00 
-1.46E+00 


Buckling  Mode  Shape  5  Critical  Load  Factor  0.4304E+02 


Mode  Shape  Contours 

X  Translation 

1.64E+01 
[  1.52E+01 

!  1.40E+01 

!  1.28E+01 

1.16E+01 
1.04E+01 
9.27E+00 
8.09E+00 
6.91  E+00 
!  5.73E+00 

4.55E+00 
!  3.37E+00 

[  2.19E+00 

j  1.01  E+00 
“  -1.75E-01 
J  -1.36E+00 
J  -2.54E+00 
J  -3.72E+00 
-4.90E+00 


Figure  4.4:  Buckling  Mode  Shapes  for  Mode  2  and  5  Predicted  by  the  20-Noded  Solid 

Elements 
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Mode  Shape  Contours 


Buckling  Mode  Shape  2  Critical  Load  Factor  0.1332E+02 


X  Translation 

3.64E+00 

3.31  E+00 
2.98E+00 
2.64E+00 

2.31  E+00 
1.98E+00 
1.65E+00 
1.32E+00 
9.85E-01 
6.53E-01 
3.21  E-01 
-1.10E-02 
-3.43E-01 
-6.75E-01 
-1.01  E+00 
-1.34E+00 
-1.67E+00 
-2.00E+00 
-2.33E+00 


Buckling  Mode  Shape  5  Critical  Load  Factor  0.3909E+02 


Mode  Shape  Contours 

X  Translation 

1.39E+01 
[  1.28E+01 

[  1.18E+01 

[  1.08E+01 

!  9.79E+00 

1  8.77E+00 

7.75E+00 
,  6.74E+00 

.  5.72E+00 

4.70E+00 
!  3.68E+00 

!  2.66E+00 

!  1.65E+00 

[  6.28E-01 

J  -3.90E-01 
J  -1.41  E+00 
J  -2.43E+00 
J  -3.44E+00 
-4.46E+00 


Figure  4.5:  Buckling  Mode  Shapes  for  Mode  2  and  5  Predicted  by  NURBS-based  Solid 
Element  Model  of  Order  3x3x2  Using  3x3x2  Numerical  Integration 
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Mode  Shape  Contours 


Buckling  Mode  Shape  2  Critical  Load  Factor  0.1814E+02 


X  Translation 


6.29E+00 

5.94E+00 

5.59E+00 

5.24E+00 

4.89E+00 

4.54E+00 

4.19E+00 

3.83E+00 

3.48E+00 

3.13E+00 

2.78E+00 

2.43E+00 

2.08E+00 

1.73E+00 

1.38E+00 

1.03E+00 

6.80E-01 

3.29E-01 

-2.13E-02 


Buckling  Mode  Shape  5  Critical  Load  Factor  0.4596E+02 


Mode  Shape  Contours 

X  Translation 

3.92E+00 
!  3.13E+00 

!  2.33E+00 

1.53E+00 
7.37E-01 
-5.95E-02 
-8.56E-01 
-1.65E+00 
-2.45E+00 
!  -3.25E+00 
-4.04E+00 
-4.84E+00 
J  -5.64E+00 
J  -6.43E+00 
!  -7.23E+00 
J  -8.03E+00 
J  -8.82E+00 
J  -9.62E+00 
-1  04E+01 


Figure  4.6:  Buckling  Mode  Shapes  for  Mode  2  and  5  Predicted  by  NURBS-based  Solid 
Element  Model  of  Order  5x5x2  Using  4x4x3  Numerical  Integration 
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Figure  4.7:  Problem  Definition  for  Buckling  of  a  Square  Plate  Subjected  to  Biaxial 
Compression  and  the  Control  Net  of  a  NURBS-based  Solid  Element  Model  of  Order 
5x5x2  Representing  One  Quarter  of  the  Plate  (Shaded  Area) 
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Mode  Shape  Contours 


Buckling  Mode  Shape  4  Critical  Load  Factor  0.5259E+09 


Z  Translation 

6.00E-01 
5.33E-01 
4.66E-01 
4.00E-01 
3.33E-01 
2.66E-01 
1.99E-01 
1.33E-01 
6.58E-02 
-9.99E-04 
-6.78E-02 
-1.35E-01 
-2.01  E-01 
-2.68E-01 
-3.35E-01 
-4.02E-01 
-4.68E-01 
-5.35E-01 
-6.02E-01 


Buckling  Mode  Shape  9  Critical  Load  Factor  0.1439E+10 


Mode  Shape  Contours 


Z  Translation 


3.66E-01 
3.25E-01 
2.85E-01 
2.45E-01 
2.04E-01 
1.64E-01 
1.24E-01 
8.32E-02 
4.28E-02 
2.43E-03 
-3.80E-02 
-7.83E-02 
-1.19E-01 
-1.59E-01 
-1.99E-01 
-2.40E-01 
-2.80E-01 
-3.21  E-01 
-3.61  E-01 


Figure  4.8:  Buckling  Modes  4  and  9  of  the  Bi-axially  Loaded  Plate  Predicted  by  NURBS- 
based  Solid  Element  Model  of  Order  5x5x2  Using  8x8x3  Numerical  Integration 
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Buckling  Load  Factor 


Figure  4.9:  Buckling  Load  Factors  Obtained  by  NURBS-based  Solid  Elements  of  Various 
Polynomial  Orders  Using  Proposed  Default  Orders  of  Numerical  Integration 
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Mode  Shape  Contours 


Buckling  Mode  Shape  2  Critical  Load  Factor  0.1 1 35E+02 


X  Translation 


2.53E+00 
2.1  6E+00 
1.80E+00 
1.43E+00 
1.07E+00 
7.00E-01 
3.35E-01 
-3.1  0E-02 
-3.97E-01 
-7.62E-01 
-1.13E+00 
-1.49E+00 
-1.86E+00 
-2.22E+00 
-2.59E+00 
-2.96E+00 
-3.32E+00 
-3.69E+00 
-4.05E+00 


Buckling  Mode  Shape  5  Critical  Load  Factor  0.3072E+02 


Mode  Shape  Contours 


X  Translation 

8.84E+00 
8.14E+00 
7.45E  +  00 
6.75E  +  00 
6.05E+00 
5.36E  +  00 
4.66E+00 
3.96E+00 
3.27E  +  00 
2.57E+00 
1.87E+00 
1.18E  +  00 
4.79E-01 
-2.17E-01 
-9.1  4E-01 
-1.61E+00 
-2.31  E+00 
-3.00E  +  00 
-3.70E+00 


Figure  4.10:  Buckling  Modes  2  and  5  of  the  Cylindrical  Shell  Predicted  by  NURBS-based 
Solid  Element  Model  of  Order  5x5x2  Using  5x5x2  Numerical  Integration 
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Displacement  Contours 


Y  Z 

Displacements;  Load  Case  1 : 


X  Translation 

2.42E-03 
-1.50E-01 
-3.03E-01 
-4.55E-01 
-6.08E-01 
-7.60E-01 
-9.13E-01 
-1.07E+00 
-1 .22E  +  00 
-1.37E  +  00 
-1.52E+00 
-1.68E+00 
-1.83E  +  00 
-1.98E+00 
-2.13E+00 
-2.29E+00 
-2.44E  +  00 
-2.59E+00 
-2.74E+00 


Contour  Plot 
Stress  Components 
Load  Case  1 
20N/Solid  SVM 


6.98E-03 
6.61  E-03 
6.25E-03 
5.88E-03 
5.51  E-03 
5.15E-03 
4.78E-03 
4.41  E-03 
4.05E-03 
3.68E-03 
3.31  E-03 
2.95E-03 
2.58E-03 
2.21  E-03 
1.85E-03 
1.48E-03 
1.1  IE-03 
7.46E-04 
3.79E-04 


Figure  4.11:  Linear  Static  Solutions  Obtained  Using  NURBS-Based  Solid  Element  Model 
of  Polynomial  Order  5x5x2  Using  4x4x3  Numerical  Integration 
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Displacement  Contours 


Y  Z 

Displacements;  LoadCase  1: 


X  Translation 

1.67E-03 
-1.53E-01 
-3.08E-01 
-4.62E-01 
-6.17E-01 
-7.72E-01 
-9.26E-01 
-1.08E  +  00 
-1.24E+00 
-1.39E  +  00 
-1.54E  +  00 
-1.70E+00 
-1.85E+00 
-2.01  E  +  00 
-2.16E  +  00 
-2.32E+00 
-2.47E  +  00 
-2.63E  +  00 
-2.78E+00 


Contour  Plot 
Stress  Components 
Load  Case  1 
20N/Solid  SVM 


7.58E-03 

7.18E-03 

6.78E-03 

6.38E-03 

5.98E-03 

5.58E-03 

5.18E-03 

4.78E-03 

4.38E-03 

3.98E-03 

3.58E-03 

3.18E-03 

2.78E-03 

2.39E-03 

1.99E-03 

1.59E-03 

1.19E-03 

7.88E-04 

3.89E-04 


Figure  4.12:  Linear  Static  Solutions  Obtained  Using  NURBS-Based  Solid  Element  Model 
of  Polynomial  Order  5x5x2  Using  5x5x2  Numerical  Integration 
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5.0  COMPARISON  OF  COMPUTATIONAL  EFFICIENCY  OF  THE  NURBS- 

BASED  AND  TRADITIONAL  SOLID  ELEMENTS  FOR  SHIP  STRUCTURE 
ANALYSES 


In  order  to  benchmark  the  computational  efficiency  of  the  NURBS-based  solid  element  against 
the  traditional  20-noded  isoparametric  solid  element  for  analysis  of  practical  ship  structures,  a 
simplified  ship  hull  component  was  considered.  This  simplified  model  included  hull  plating,  a 
truncated  transverse  bulkhead  and  two  stiffeners  with  Tee  section,  as  shown  in  Figure  5.1.  The 
original  model  was  a  surface  model.  In  the  present  work,  this  surface  model  was  first  reverse- 
engineered  utilizing  Intelliship  (Intergraph  Corporation,  www.intergraph.com)  to  create  a  solid 
model. 

The  solid  model  was  created  with  help  of  Rhino3D  (McNeel  Corporation,  www.rhino3d.com). 
The  procedures  used  to  create  a  solid  model  based  on  the  surface  model  include  the  following 
steps: 

1 .  Offset  each  surface  in  the  surface  model  by  the  distances  determined  by  the  structural 
object  thickness  to  generate  surfaces  representing  top,  bottom  and  interior  of  the 
structural  components. 

2.  Trim  and  extend  the  offset  surfaces  at  the  intersections  of  structural  components. 

3.  Create  new  surfaces  for  the  solid  entities  caused  by  the  thicknesses  of  the  intersecting 
components. 

4.  Sew  all  the  surfaces  to  form  solids. 

After  the  solid  model  was  completed,  it  was  used  to  create  an  isogeometric  NURBS  model, 
which  can  be  defined  as: 

•  It  must  be  a  solid  model  with  6  surfaces. 

•  Two  opposite  surfaces  must  have  the  same  NURBS  definitions,  e.g.  top  and  bottom 
surfaces  must  have  the  same  degree,  number  of  control  points  and  same  knot  vectors. 

•  If  two  solid  patches  are  connected,  then  the  two  surfaces  that  touch  each  other  must 
have  identical  NURBS  definitions. 

In  order  to  create  a  NURBS  model,  which  satisfies  all  the  above  requirements,  one  command 
for  checking  the  surface  parametric  directions  is  implemented  for  Rhino3d.  During  the  model 
generation,  a  new  set  of  surface  definitions  was  recreated  for  all  the  surfaces  in  the  solid  model 
so  that  all  the  requirements  for  NURBS  model  definitions  were  fulfilled.  Then  with  help  of  the 
Rhino3d  command,  the  orientation  of  the  surface  parametric  directions  was  fixed.  The  final 
solid  model  generated  using  this  approach  is  shown  in  Figure  5.2.  This  model  includes 
NURBS  representation  of  each  individual  solid  patch  involved  in  the  finite  element  analysis. 
However,  to  obtain  a  NURBS-based  solid  element  model  useable  for  stress  analysis,  node 
equivalencing  must  still  be  performed  to  properly  connect  these  solid  patches.  A  special 
purpose  FORTRAN  program  was  developed  to  perform  this  task  of  model  integration. 

There  were  some  issues  encountered  during  the  creation  of  the  solid  model.  The  most 
important  ones  are  presented  below: 
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•  Rhino3d  gives  errors  when  trying  to  trim  the  hull  with  bulkhead  surfaces  (or  even  with 
a  plane).  This  problem  is  avoided  by  increasing  the  thickness  from  7mm  to  10mm.  This 
is  a  very  serious  problem  because  the  CAD  tool  cannot  handle  the  required  geometry. 
One  solution  is  to  implement  a  new  application  to  do  the  trimming  using  third  party 
geometric  kernel  software,  such  as  TSNLib  (Spatial  Corporation)  or  ACIS  (SML 
Corp),  export  the  model  into  standard  file  formats  like  IGES,  and  then  import  the 
model  into  the  CAD  tool.  However,  this  is  a  time  consuming  procedure. 

•  When  creating  a  NURBS  representation  of  a  solid  entity,  any  of  the  commercial  CAD 
tools  provide  user  controls  for  the  math  definitions  of  the  surfaces  that  define  a  solid 
patch.  These  math  definitions  of  surfaces  include  polynomial  orders,  knot  vectors  and 
number  and  distribution  of  the  control  points.  However,  in  a  NURBS-based  solid 
element  model,  not  only  must  the  math  definitions  of  opposite  surfaces  in  a  given  solid 
patch  be  identical,  but  the  math  definitions  of  the  contiguous  surfaces  between  two 
connecting  solid  patches  must  also  be  identical  as  well.  In  order  to  resolve  this 
difficulty,  separate  tools  were  developed  in  the  present  work.  However,  the  tools  that 
were  developed  for  this  work  are  neither  general-purpose  nor  based  on  any  commercial 
CAD  tools.  To  create  models  for  the  NURBS-based  solid  element  efficiently,  an 
automatic  tool  based  on  commercial  CAD  programs  like  Rhino3d  should  be  developed. 

To  perfonn  linear  static  analysis,  clamped  boundary  conditions  were  applied  to  the  two  sides 
and  the  bottom  edge  of  the  model.  A  concentrated  force  were  applied  close  to  the  intersection 
of  the  bulkhead  and  one  of  the  stiffeners  on  the  outside  of  the  ship  hull  pointing  inwards  in  the 
direction  normal  to  the  surface  as  depicted  in  Figure  5.3.  Although  these  selections  of  the 
boundary  and  loading  conditions  did  not  correspond  to  a  practical  loading  case  on  ship 
structures,  they  provided  sufficient  constraints  and  loads  to  the  numerical  model  to  pennit 
comparison  of  the  local  deformations  and  stresses  at  the  junctions  of  different  structural 
components  predicted  by  various  finite  element  formulations.  The  undeformed  control  net  and 
the  deformed  configuration  of  the  stiffened  panel  are  shown  in  Figure  5.4,  where  a  summary  of 
the  displacements  at  one  of  the  loaded  nodes  is  also  shown.  This  plot  of  deformed  shape  was 
generated  using  the  post-processing  capability  for  the  NURBS-based  finite  element  developed 
in  this  contract. 

In  order  to  compare  the  computational  efficiency  of  the  NURBS-based  finite  element  and 
traditional  20-noded  isoparametric  solid  element,  two  meshes  of  the  20-noded  solid  elements 
were  constructed  for  the  same  stiffened  plate  component  with  the  aid  of  the  post-processing 
capability  for  NURBS-based  element.  In  one  of  the  solid  element  meshes,  each  element  in  the 
NURBS-based  solid  element  model,  which  is  defined  by  knot  intervals  in  the  parametric 
space,  was  converted  into  one  20-noded  solid  element.  This  mesh  is  shown  in  Figure  5.5.  In 
the  second  solid  element  mesh,  each  element  in  the  NURBS-based  model  was  subdivided  into 
eight  20-noded  solid  elements  as  shown  in  Figure  5.6.  Linear  static  analyses  were  performed 
using  both  solid  element  models.  The  deformed  configurations  predicted  by  the  solid  elements 
are  given  in  Figures  5.5  and  5.6  along  with  summaries  of  nodal  displacements.  A  comparison 
of  the  displacements  indicated  that  the  accuracy  of  the  solutions  from  the  NURBS-based  solid 
element  model  is  comparable  to  those  obtained  using  the  finer  mesh  of  the  20-noded  solid 
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element.  However,  as  indicated  in  Table  5.1,  the  computation  time  required  by  the  NURBS- 
based  model  is  even  less  than  that  required  by  the  coarse  mesh  of  the  20-noded  solid  element. 

In  order  to  further  evaluate  the  relative  performance  of  the  NURBS-based  and  20-noded  solid 
elements,  the  von  Mises  stress  contours  generated  by  the  various  models  were  compared.  The 
stress  distribution  predicted  by  the  NURBS-based  elements  is  shown  in  Figure  5.7,  where  both 
the  entire  model  and  the  local  details  in  the  vicinity  of  the  applied  loads  are  displayed.  The 
solutions  generated  by  the  two  20-noded  solid  element  meshes  for  the  same  local  area  are 
presented  in  Figure  5.8.  A  comparison  of  these  results  revealed  that  the  stress  solution  from  the 
NURBS-based  solid  element  model  has  a  significantly  higher  resolution  than  the  results  of  the 
coarse  solid  element  model. 

These  results  clearly  demonstrated  the  potential  of  the  NURBS-based  finite  element  for  ship 
structure  analyses. 
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Table  5.1:  Comparison  of  Computation  Times  Required  by  the  NURBS-based  and  20- 

Noded  Solid  Finite  Element  Models 


Model 

No.  of 
Nodes 

Elapsed  Time  By  Mod 

ule  (sec) 

Total  Time 
(sec) 

Ratio 

ELEMS1 

DECOM 

STRESS 

NURBS-based 

5532 

14.297 

12.218 

1.594 

28.594 

1:1 

Solid  (coarse) 

12071 

7.953 

30.500 

0.656 

39.453 

1.38:1 

Solid  (fine) 

74513 

41.594 

774.593 

41.859 

968.140 

33.66:1 
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Figure  5.1:  Surface  Model  of  a  Stiffened  Plate  Structural  Component 


Figure  5.2:  Isogeometric  NURBS  Model  of  the  Stiffened  Plate  Structural  Component 
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Figure  5.3:  Control  Net  of  the  NURBS-based  Solid  Element  Model  of  the  Stiffened  Plate 
Structural  Component  Shown  with  Boundary  Conditions  and  External  Load 
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Figure  5.4:  Undeformed  Control  Net  and  Deformed  Configuration  of  a  NURBS-based 
Solid  Element  Model  of  the  Stiffened  Plate  Structural  Component 
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Figure  5.5:  Undeformed  and  Deformed  Meshes  of  the  Coarse  Mesh  of  20-Noded  Solid 
Element  of  the  Stiffened  Plate  Structural  Component 
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Figure  5.6:  Undeformed  and  Deformed  Meshes  of  the  Fine  Mesh  of  20-Noded  Solid 
Element  of  the  Stiffened  Plate  Structural  Component 
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Figure  5.7:  Stress  Contour  Predicted  by  the  NURBS-based  Solid  Element  Model  of  the 

Stiffened  Plate  Structural  Component 
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Figure  5.8:  Stress  Contour  Predicted  by  the  Coarse  (Top)  and  Fine  (Bottom)  Solid 
Element  Models  of  the  Stiffened  Plate  Structural  Component 
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CONCLUSIONS 


In  this  report,  we  presented  work  perfonned  in  the  second  phase  of  a  project  for  developing  a 
NURBS-based  computational  method  for  the  VAST  finite  element  program.  Within  the  initial 
phase  of  the  development,  a  NURBS-based  solid  element  was  implemented  and  extensively 
verified  for  various  types  of  linear  and  nonlinear  analyses  [1],  In  the  present  phase,  the  earlier 
development  was  continued  in  a  number  of  directions,  including  refinement  of  the 
convergence  study,  investigation  of  the  accuracy  problems  for  buckling  analyses,  evaluation  of 
computational  efficiency  and  development  of  a  post-processing  capability.  All  these  objectives 
were  completed  successfully. 

In  the  convergence  study  performed  in  the  previous  phase,  the  reference  solution  was  obtained 
using  a  relatively  coarse  mesh  of  20-noded  solid  elements,  so  the  accuracy  of  the  reference 
solution  itself  was  questionable.  In  the  present  study,  additional  results  for  20-noded  solid 
elements  were  obtained  using  various  mesh  layouts  and  numerical  integration  orders.  This 
extended  convergence  study  confirmed  that  the  NURBS-based  solid  element  converges  to  the 
correct  solutions  with  a  convergence  rate  that  is  much  faster  than  that  of  the  traditional  20- 
noded  isoparametric  solid  element. 

The  investigation  on  the  accuracy  of  buckling  analysis  using  the  NURBS-based  solid  element 
was  started  from  consideration  of  a  simple  test  case  involving  a  simply-supported  square  plate 
subjected  to  uniform  bi-axial  compressive  stress.  This  test  problem  was  selected  because  of  the 
availability  of  analytical  solutions.  Numerical  solutions  to  this  problem  were  obtained  using 
various  NURBS-based  models  and  different  orders  of  numerical  integrations.  The  numerical 
results  were  then  compared  with  the  analytical  solution.  These  comparisons  indicated  that  in 
order  to  obtain  reliable  and  accurate  results  for  linear  buckling  analyses,  the  orders  of  the 
numerical  integration  must  be  at  least  equal  to  the  polynomial  orders  of  the  NURBS  basis 
functions.  This  has  been  implemented  as  the  default  integration  orders  in  the  VAST  program. 
Using  this  default  rule  for  numerical  integration,  the  test  problem  considered  in  the  previous 
phase  was  reconsidered  and  very  good  convergence  properties  were  obtained. 

In  order  to  establish  a  benchmark  for  computational  efficiency  of  the  NURBS-based  element 
compared  with  the  traditional  20-noded  isoparametric  solid  element,  a  simplified  structural 
model  consisting  of  a  hull  plate  and  intersecting  truncated  bulkhead  and  two  longitudinal 
stiffeners  was  considered.  The  construction  of  the  NURBS-based  model  required  reverse¬ 
engineering  of  the  CAD  surface  model.  This  reverse-engineering  involved  the  following 
operations:  offsetting  and  duplicating  surfaces,  modifying  surface  definitions  to  ensure  that  the 
intersections  between  solid  entities  are  generated  properly  and  generating  fully  compatible 
NURBS  representations  at  interfaces  of  opposite  surfaces.  A  special-purpose  computer  tool 
was  developed  for  this  task.  The  NURBS-based  model  of  the  stiffened  plate  structural 
component  was  used  in  a  linear  static  analysis  and  the  results  were  compared  with  two  20- 
noded  solid  element  models  with  different  levels  of  mesh  refinement.  The  benchmark  results 
indicated  that  in  order  to  achieve  comparable  accuracy,  the  NURBS-based  element  requires  a 
much  coarser  mesh,  leading  to  a  significant  saving  of  computation  time. 
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The  results  from  the  present  study  further  confirmed  the  potential  of  the  NURBS-based  solid 
element  in  engineering  analyses  and  also  the  study  significantly  improved  its  readiness  for 
practical  applications. 
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